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Abstract. Lattice four-fermion models containing N flavors of staggered fermions, that are invariant 
under Z2 and 17(1) chiral symmetries, are known to suffer from sign problems when formulated using 
the auxiliary field approach. Although these problems have been ignored in previous studies, they can 
be severe. In this talk, we show that the sign problems disappear when the models are formulated in the 
fermion bag approach, allowing us to solve them rigorously for the first time. 



1. Introduction 

Four-fermion field theories are interesting in both condensed matter and particle physics. The well known 
Hubbard model and its variants are often used in studying cuprate superconductors [ 1], antiferromagnets 
[2] and more recently graphene [3]. Low energy nuclear physics is also studied with four-fermion 
couplings in the effective field theory framework (4j Q . In the context of more fundamental theories 
like QCD, four-fermion field theories offer a simpler setting to study phenomena like fermion mass 
generation and chiral symmetry breaking It has been suggested recently that quantum critical 
phenomena in graphene can be studied with four-fermion field theories Q[8]]. It has also been found 
there is a QCD-like sign problem in the four-fermion field theory Despite the wide interest, strongly 
coupled four-fermion field theories remain poorly understood as compared to their bosonic counterparts 
due to computational difficulties. 

The only available method to compute quantities in a strongly interacting field theory with no small 
parameter is the Monte Carlo (MC) method. Due to the quantum nature of a fermion there are no 
natural fermion configurations with positive weights that can be used for important sampling. In two 
space-time dimensions fermions can often be bosonized and models can we written in terms of world 
line configurations with positive weights. This fact can be used to design powerful MC methods 
ifTUl [TT1 [T2l . In higher dimensions, the traditional MC approach is to integrate the fermions out in 
favor of a determinant of a large fermion matrix. Whenever this determinant is positive a non-local 
probability distribution emerges, which can be used to construct a MC method. The most popular is 
the Hybrid Monte Carlo (HMC) method lfT3l [141 which has continued to evolve in many ways since 
its discovery llT5l . Unfortunately, small eigenvalues of the fermion matrix which naturally arise in the 
presence of massless fermions can cause singularities in the HMC approach. This makes it difficult 
to study quantum critical phenomena containing massless fermions. While other determinantal MC 
methods do not encounter such problems, they scale poorly with system size iPToTl . In cases where the 
determinant of the fermion matrix is not positive, the original theory is said to suffer from a sign problem 
and the traditional approach is not useful. The repulsive Hubbard model away from half filling is a 
classic example where progress has been limited due to sign problems. Other relativistic four-fermion 



field theories like the Gross-Neveu (GN) models and Nambu-Jona-Lasinio (NJL) models are also known 
to suffer from sign problems in three or more space-time dimensions ifTTl . 

Recently a new approach called the fermion bag approach was proposed to solve some four-fermion 
field theories |fT8l[T9ll20ll2Tll22ll . It is an extension of the meron cluster idea proposed some time ago 
If23l . The idea behind the fermion bag is to identify fermion degrees of freedom that cause sign problems 
and collect them in a bag and sum only over them. This is in contrast to traditional approaches where all 
fermion degrees of freedom in the entire thermodynamic volume are summed to solve the sign problem. 
When the fermion bag contains only a small fraction of all the degrees of freedom and the summation 
can be performed quickly, the fermion bag approach can be used to design powerful MC methods. 
Sometimes, the bag splits into many disconnected pieces further simplifying the calculation. The fermion 
bag approach has three main advantages: (a) Due to a duality, fermion bag sizes are small both at weak 
and strong couplings, (b) Singularities in the massless limit can be tackled without a problem, (c) Some 
sign problems that haunt traditional approaches are naturally solved. While the first two advantages 
have been demonstrated, the third advantage is not so clear from previous work. Here we show how 
solutions to some unsolved sign problems in four-fermion models also emerge naturally in the fermion 
bag approach. 

It is useful to clarify some confusions that may arise about what we mean by a sign problem and thus 
a solution to the sign problem. If one can write the partition of a quantum statistical mechanics system 
as a sum over configurations whose Boltzmann weights are all positive and if the cost of computation 
of the Boltzmann weights only scales as a polynomial in system size, then we say the model does not 
suffer from a sign problem. However, as already stated above, in fermionic systems there are no natural 
configurations where the Boltzmann weights are positive. The conventional method is to use the auxiliary 
field approach to expand the partition function as a sum of bosonic configurations where the fermion 
determinant is taken as part of the Boltzmann weight. If this weight can be negative one often says 
the model suffers from a sign problem. However, if an alternate approach can be found where the sign 
problem disappears, one can of course say the problem never suffered from a sign problem to begin 
with. On the other hand, if this alternate approach was not known earlier, the new approach can be 
considered as a solution to the sign problem present in the other method. This is what we mean when 
we say "solutions to unsolved sign problems". It must be noted that all sign problems are problems in 
exactly this sense. Once a solution is found there is no longer a problem. It is of course likely that some 
problems may remain unsolved j24l . 

We consider lattice GN models containing N flavors of massless staggered fermions with either a Z2 
or a £7(1) chiral symmetry IfTTl . While we work in three space-time dimensions, our results can easily 
be extended to higher dimensions. Although in three dimensions the symmetries we refer to are a part 
of a flavor symmetry, they are often loosely called chiral symmetries in the literature. The Z2 models 
with odd N and all the £7(1) models are known to suffer from a sign problem when formulated in the 
traditional auxiliary field approach. Here we show that the sign problems disappear in the fermion bag 
approach. Our paper is organized as follows. In section 2 we review the auxiliary field approach to 
lattice GN models with both Z2 and £7(1) chiral symmetries and discuss how the sign problems arise. In 
section 3 we discuss the severity of the sign problems. In section 4 we discuss the fermion bag approach 
and show that sign problems do not arise. Section 5 contains our conclusions. 

2. Auxiliary field approach 

Lattice GN models are formulated in the auxiliary field approach through the action 



where Xi{ x )i Xi{ x ) denote the Grassmann valued fermion fields of flavor i = 1, 2.., N at the lattice site 
x. The explicit form of the auxiliary field action Saf depends on the GN model and will be discussed 




(1) 



below. The matrix D [(f)] is defined by 

(D{4>]) xy = D xy + 5 xy 0(x), (2) 

where 4>(x) is a function of the auxiliary fields as defined below and D x>y is the free staggered fermion 
matrix 11231126112711. 

o=l,2,3 

Since we work in three dimensions, a labels the three directions, r\ x ^ a = e^^ a ' x \ (1 = (0,0,0), 
£2 = (1, 0, 0), £3 = (1,1,0) are the staggered fermion phase factors and m is the bare fermion mass. We 
assume anti-periodic boundary conditions in all directions and denote the lattice volume by V = L i . 

Following Cm . we define the auxiliary fields on dual sites x. The model with a Z2 chiral symmetry 
is defined through a single real auxiliary field a(x), such that 



" x 

(x,x) 

while the model with a U(l) chiral symmetry requires two real auxiliary fields a(x) and tt(x), such 
that 
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where e(x) is the parity of a lattice site (1 on even sites and —1 on odd sites). In the above expressions, 
the set of nearest dual sites x, surrounding the fixed lattice site x is denoted as (x, x) (see Fig. [I]). In this 
work we only consider these two classes of models. 
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Figure 1. Nearest neighbor lattice sites (open circles) of a fixed dual site (filled circle) x is represented 
by [x, x] (see left figure), while the nearest neighbor dual sites of a fixed lattice site is denoted by (x, x) 
(see right figure). 

It is easy to verify that Sqn is invariant under U(N) flavor transformations. When m = 0, 
additional chiral symmetries emerge. The Z2 model is invariant under Xi( x ) ~^ £ ( x )Xi( x )i Xi( x ) ~^ 
—Xi( x )s( x ), cr(i) — > — &( x ) while the U(l) model is invariant under the additional U(l) chiral 



20 40 60 80 100120140 
log|Det(D[0])| 



20 40 60 80 100120140 
log|Det(D[0])| 




550 600 650 700 

log]Det(D$])] 




550 600 650 700 
log]Det(D[?l)] 



Figure 2. Distributions of positive (left graphs) and negative (right graphs) weight configurations as a 
function of log |Det(D[0])|. One million configurations and 5000 configurations were generated at 6 3 
and 12 3 lattices respectively. The distribution of positive configurations is almost identical to the one 
with negative configurations suggesting a severe sign problem. 

symmetry Xi(x) e* e ^^/ 2 Xi(^), Xii x ) — > Xiix)^^ 9 ^ 2 , cr(x) ^ a (x) cos 9 + ir(x) sin 6, ir(x) — > 
7r(x)cos6> — a(x)sin9. The models contain a quantum critical point (QCP) separating a chirally 
symmetric phase (at small couplings) from a phase where the chiral symmetry is spontaneously broken 
(at large couplings). The symmetries that govern the QCP needs proper analysis due to fermion doubling. 
Without such an analysis it is difficult to establish the continuum field theory that emerges at the critical 
point EOT 

In the traditional MC approach, one integrates over the Grassmann fields and writes the partition 
function of the GN models as 
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[VaVn] e~ SAF[a ^ ] \ DetD([0]) 



N 



(8) 
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In order to design a MC method the determinant terms in the above expressions have to be real and 
positive. In the Z2 model since is real, the matrix elements of D[4>] are real. Hence, the determinant is 
real but not necessarily positive. In the case of the U(l) model, 4> is complex and so the matrix elements 
of D [4>] and its determinant can be complex. Hence, the Z2 model as formulated in Eq. ( 8 1 suffers from 
a sign problem for all odd values of N, while the model as formulated through Eq. (9 1 suffers from 
a sign problem for all values of N. 



3. Severity of the sign problem 

Earlier calculations in the Z2 and U(l) GN models have all been performed in the auxiliary field 
approach |[28l l29l [30l . The essential focus has been to understand the quantum phase transition and 
compute the critical exponents. These calculations have circumvented the sign problem by studying 
even N in the Z2 case or by introducing conjugate fermions with an opposite chiral charge in the £7(1). 



Inclusion of conjugate fermions changes the partition function from Eq. Q to 



Z c u{ i) = /MM e~ SAF 



DetD([0]) 



2jV 



(10) 



and changes the flavor symmetries to U (JV) x U(N) while the chiral symmetry remains unchanged. 

The TV = 1 model with Z2 chiral symmetry was also studied in the auxiliary field approach using 
the HMC algorithm (31). Strangely, in this study the sign of the determinant was never discussed and 
seems to have been ignored. Since the results of the quantum critical behavior were in quantitatively 
agreement with large N results (improved with Pade-approximations), it may have been assumed that 
the sign problem was mild. If this is indeed true then statistically, positive sign configurations should 
dominate over negative sign configurations. The Z2 model studied in iPTI is slightly different from the 
model studied here. The auxiliary fields a(x) also live on the main lattice site and the field (ft appearing 
in the Dirac operator D[<j>] of Eq. (|2j), is defined as (ft{x) = | Yl(z x) a ( z )' where now { z i %) refers to the 
six nearest neighbor sites z for a given site x. In order to study the sign problem, we generated several 
Gaussian random auxiliary field configurations according to the distribution 

P(a(x)) = exp ^ - " I log(vr)}^ (11) 

and computed Det(D[4>]) for each of these configurations. We then separated the configurations into 
those with a positive determinant and those with a negative determinant. In Fig. [2] we plot the distribution 
of configurations with positive and negative determinants as a function of log |Det(D[0])| for 6 3 and 12 3 
lattices. As can be seen, the distribution of configurations with positive and negative weights are almost 
identical suggesting a severe sign problem rather than a mild one! Although we are not performing 
important sampling, our results clearly show that the sign problem must be studied carefully. 

An important question to study is whether the HMC algorithm is getting trapped in the sector of 
configurations with positive weights (or negative weights). Note that, in the Z2 models the only way to 
move from a positive weight sector to the negative weight sector is to pass through configurations which 
have almost zero weight assuming the step size in the HMC algorithm is small. Perhaps the suppression 
of the tunneling between the two sectors leads to long auto-correlation times or even lack of ergodicity. 
This argument also applies to Z2 models with even 



4. Fermion bag approach 

We will now show that the sign problems in both the Z2 and the U(l) models discussed in section [2] 
disappear in the fermion bag approach. The proof relies on the fact that any k{ -point correlation function 
involving the i th flavor of staggered fermions defined through 

Ci{x h ,...,x ik ,) = J [dXid X i]e~^y * i{x) D *» Xi{y) Xi{xiMx h ) ...Xi(xi ki )Xi(xi ki ) (12) 

is positive semi-definite. This is due to the special properties of the free staggered fermion matrix. 
Indeed, using the ideas developed in the fermion bag approach EOl . we can write 

Ci( Xil ,..,x ik J = Det(D) Det(G[{x}i]) = Det(W [{*}*]) (13) 

where is the ki x k% matrix of propagators between the k% sites in the set {x}i = Xi p , p = 1, .., ki 

whose matrix elements are G Xp)Xq = D~ l x and the matrix Vl^[{x}i] is a (V — ki) x (V — ki) matrix 
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identical to the matrix D except that the sites in the set {x}i are dropped from the matrix. All the 



determinants appearing in Eq.( 13 ) can be shown to be positive (or zero). The simplest way to see this is 
to consider the matrix W. Since it is exactly the same as the staggered fermion matrix with some sites 
removed, its eigenvalues come in complex conjugate pairs of the form m ± iX. Unpaired eigenvalues 
are always m and they too come in pairs when the lattice is bipartite. When m = then the determinant 
can be exactly zero. Thus, Ci(xi 1 , .., Xi k . ) > 0. We will use this property to prove the absence of a sign 
problem in the fermion bag approach. 

Instead of integrating out the fermion fields let us integrate out the auxiliary fields first and construct 
the appropriate four fermion action for the models. Let us first consider the Z2 model. Each integral over 
the auxiliary field a(x) on the dual site x gives, 



(14) 



where M = y / 2irg 2 /N and 



Si(x) 
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128N 



(15) 



is the effective four-fermion interaction term at each dual site x. The symbol [x, x] denotes the set of all 
lattice sites surrounding the dual site x (see Fig. [T]) . Thus, each integral generates many four-fermion 
couplings of the form Xi( x )Xi( x )Xj(y)Xj(y) where i and j are arbitrary flavor indices and x and y are 
corners of the cube surrounding the dual site x. We can classify the possible couplings into four types 
based on the bonds (xy) connecting the corners x and y. If the two corners are the same we refer to 
it as a site-bond or a 5-bond. If the two corners are the two neighboring sites we get a L-bond (or a 
link-bond). Similarly, if the two corners are across a face diagonal or a body diagonal, we call the bonds 
F-bond and S-bond respectively. These four bond types are illustrated Fig. [3] 




Figure 3. An illustration of the four types of four-fermion couplings (or bonds) generated through the 
auxiliary field integration. From left to right we have a S, L, F and B bond respectively. 

Integration over all the auxiliary field variables yields the four-fermion interaction term of the action 
Sz 2 ,mt = Ylx Collecting the terms in each of the four types of four fermion couplings separately 

we see that 

Sz 2 ,int = U S B S + U L B L + U F B F + U B B B (16) 
where C/ s /4 = U L /4 = U F /2 = U B = 5 2 /(64iV) and 

^bond = Xi( x )Xi(x)Xj(y)Xj(y)- (17) 

i,j',(a;j/}Sbond 

Based on the above results, the partition function of the Z2 model can be rewritten as 

Zz 2 = J H[dXidXi]e- Sz *. (18) 



where Sz 2 = So + Sz 2 , mt is the equivalent four-fermion action of the model. Here Sq = 
T, x ,y,i Xi(x)D x , y Xi{y) is the free fermion action. 

In the fermion bag approach, each four-fermion coupling is represented as a bond and expanded in 
powers of the coupling. For example the four-fermion coupling of the type Xi( x p)Xi( x p)Xj{ x q)Xj( x q) 
can be denoted by the bond variable bij(x p , x q ) = 0, 1, such that if it is then no bond is assumed to 
exist between the sites x p and x q , otherwise the specific four-fermion coupling is inserted in the partition 
function. Due to the Grassmann nature of the couplings higher powers of the couplings do not exist. 
More details can be found in |[T8l . Thus, in the fermion bag formulation, the partition function can be 
written as a sum over these bond configurations [b] , such that 



z Z2 = Y. u s s u n L L u n /u n B B [ nfe^] e ~ 5 °n^(^)^(^)--^K)^K) 

J i i 

; ur ul L up ul B { n Ci ( Xil , . . , Xlki ) } d9) 

i 

where ns,riL,nF and ub are the total number of bonds of each type and the correlation function 



[b] 



Ci(xi 1 , .., Xi k _) was defined in Eq.( 12 1. A given bond configuration [b] uniquely determines the ki sites 
Xi 1 ....Xi k , (ordered in a consistent way). Since we argued above that Cj(xj 1 , .., xi h .) > there is no sign 
problem in this expansion of the partition function for all non-negative values of Us,Ul,Uf, Ub, any 
positive integer N and real mass m. 

In the case of the U(l) model, we need to integrate over both the auxiliary fields <t{x),tt{x) on every 
dual site. It is straightforward to verify that 

[da(x)dn(x)] e-^-^to.!*."] U*)xi(*)) 
x e-^^&Ms^MxiMxiOz)) = j\f e -Si(x) (20) 

where N = (4irg 2 /N) and 

Sj(x) = [ Xi( x )Xi( x ) ~ [ £ ( x )Xi( x )Xi(z 



(21) 



l,\X,X\ 



Interestingly, the four-fermion couplings of the type S and F get canceled between the two terms in the 
above equation. On the other hand couplings of the type L and B survive so that the four-fermion action 
for the U(l) model turns out to be 

Su(i) =S + U L B L + U B B B (22) 

with Ul/4 = Ub = g 2 /(l6N). Thus, the only difference between the Z2 and U(l) models is that the 
couplings Us = Uf = in the U(l) model. Indeed these couplings break the U(l) symmetry to a Z2 
symmetry as can be easily verified. Since we already proved that the sign problem in the Z2 model was 
absent for all non-negative values of Us, Ul, Uf, Ub and N in the fermion bag formulation, the same is 
true for the U(l) model as well. 



5. Results for U(l ) model 

It has been pointed out in ll32l that the Thirring model is different from the Gross-Nevue(GN) model in 
many ways. The critical exponents in a lattice GN model with a Uf(l) x Z2 symmetry with staggered 
fermions have been computed and it was found that v = 1.00(4) and 77 = 0.754(8) OH . In a U(l) 
lattice GN model simulation 11301 , 77 = 0.904(50) has been found. However, those simulations are all 



based on the conventional Monte Carlo methods with auxiliary fields and an extra scale m in the action. 
It could be tricky to carry out the chiral extrapolation analysis while simulations are closed to U c . The 
most importantly, as we have pointed out in the previous sections, the original theory which has U (1) 
flavor symmetry suffers from the sign problem, the conventional methods instead studied the modified 
theory which has U (1) x U(l) flavor symmetry with unchanged chiral symmetry. The expanding of the 
flavor symmetry can potentially affect the critical exponents and turn the model into different universality 
class. It is desirable to study the original model using a complete new method. 

In our previous paper ll33l . we used fermion bag approach to study the quantum critical exponents of 
the Thirring model that has the action 

Srhirring = Sq + Uh^L (23) 

It has been shown that our results agree with the early studies by the conventional approaches, but our 
data are more accurate. By formulating the GN model in the fermion bag approach, in fact we can 
show that on the lattice the difference between the two models is an additional four-fermion coupling 
along the body diagonal after integrating out the auxiliary fields. It has the action in Eq. ( [22] ) with the 
extra coupling Ub = Ul/4:, and the original theory can be simulated without sign problem under the 
fermion bag approach. The simulation of £7(1) GN model becomes straight forward, we should be able 
to compare the Thirring model to GN model at the transition with fairly small amount of modifications 
of the previous simulation. Before we start our simulations, as one can see, the body diagonal coupling 
is just quarter of the link coupling, it also preserves the U(l) chiral symmetry. We expect this extra term 
should not change the universality class of the GN model compared to the Thirring model, both models 
can have the same critical exponents. In the simulations, we actually let Ul = Ub = U, we believe this 
little modification should not have the affect on the critical exponents. 

In order to study the critical exponent, we focused on three observables (Let L be the lattice size): 
The chiral condensate susceptibility 

X = ^3 ^2(-(p x ip x ip y ipy), (24) 
x,y 



the chiral winding number susceptibility 



= (25) 



, and the ratio of fermion two-point correlator 

R f = C F (L/2-l)/C F (l), (26) 



The details of those observables are explained in lTT9ll . Since the fermions are exactly massless, in 
the vicinity of U c we expect these three observables to satisfy the following simple finite size scaling 
relations: 



x -iL 2 -n = J2fk[(U-U c )L¥ 

fc=0 
3 

<9?>=J>* [(U-U C )L$ 

3 

RfL 2+v*=J2p k [(U-U e )Lh 

k=0 



(27) 
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Figure 4. Plots of \ 1 L 2 r ', and RfL 2+v ^ as a function of U for L from 16 to 40. The solid lines 
show the combined fit which give U c = 0.0909(1), f = 0.88(1),// = 0.63(1) and rj^ = 0.37(1) with 
X 2 /d.o.f = 0.89 



where we have kept the first four terms in the Taylor series of the corresponding analytic functions. Our 
goal is to compute the critical exponents 77, v and at the quantum critical point. Plots of our data is 
shown in Fig. [4] Since our data fits very well to the expected scaling form for a whole range of lattice 
sizes, we feel confident that the corrections to scaling are small. 

We compared the critical exponents from U(l) GN model to the Thirring model in Table [I] There 
are slight deviations between two models. However, by combining data sets from two models, we can 
fit the same scaling relations quite well with x 2 /d.o.f = 1.19. The joined fitting results are listed as 
"Combined" in the table. It is clear that those critical exponents are completely different from mean-field 
analysis and early MC studies by conventional methods QUI . Our data shows that U(l) GN model and 
Thirring model belong to the same universality class. It will be very interesting to carry out the same 
study on the Z2 GN model and compare to the early results PTll . Since the J7(l) chiral symmetry has 
been broken explicitly by the face diagonal bond, the critical behaviors can be different from the U(l) 
model. 
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X 2 /d.o.f. 


Thirring 


0.85(1) 


0.65(1) 


0.37(1) 


1.3 


U{1) GN 


0.88(1) 


0.63(1) 


0.37(1) 


0.89 


Combined 


0.85(1) 


0.64(1) 


0.37(1) 


1.19 



Table 1. Results from the fit of the data to Eqs. ( 27 1. In addition to fit Thirring and U(l) model separately, 



we also combined two data sets and fit all the critical exponents using the same scaling relations. 



6. Conclusions 

The fermion bag approach provides an alternative approach to fermion field theories where solutions 
to new sign problems emerge naturally. Here we have demonstrated that some sign problems in the 
auxiliary field formulation of GN models, especially with Z2 and U(l) chiral symmetries, disappear in 
the fermion bag approach. By using the fermion bag approach, we have also shown the 17(1) GN model 
and Thirring model belong to the same universality classes. While we have not shown here, we can 
solve sign problems in some lattice field theories containing both dynamical boson and fermion fields 
with similar chiral symmetries. In these more complex models, the solutions emerge when bosons are 



formulated in the world-line approach and the fermions are formulated in the bag approach. Such an 
approach to quantum field theories was proposed in |[34ll . 

Sign problems in other fermion models with more complex symmetries are also solvable in the 
fermion bag approach. However, in many interesting cases the Boltzmann weight of a fermion bag, 
although non-negative, turns out to be a fermionant instead of a determinant [35]. Since the computation 
of the fermionant can be exponentially hard, the fermion bag approach loses its practical appeal in such 
cases. Still, we believe that there are many other interesting models where the weight of the fermion bag 
continues to be positive and computable with polynomial effort. 



Acknowledgments 

I would like to thank my collaborator S. Chandrasekharan for the wonderful collaboration and fruitful 
discussions. This work was supported in part by the Department of Energy grants DE-FG02-05ER41368 
and DE-FG02-00ER41 132. 



References 



[1 

[2 
[3 
[4 
[5 
[6 
[7 
[8 
[9 
[10 
[11 
[12 
[13 
[14 
[15 
[16 

[17 

[18 
[19 
L20 

[21 
L22 
L23 
L24 
[25 
[26 
[27 
[28 
L29 
L30 
[31 
[32 
[33 

L34 
[35 



Franz M and Tesanovic Z 2001 Phys. Rev. Lett. 87(25) 257003 

Ran Y, Hermele M, Lee P A and Wen X G 2007 Phys. Rev. Lett. 98(1 1) 1 17205 

Meng Z Y, Lang T C, Wessel S, Assaad F F and Muramatsu A 2010 Nature 464 847 

Bedaque P F and van Kolck U 2002 Ann.Rev.Nucl.Part.Sci. 52 339-396 (Prepnnf |r7ucl-th/0203055) 

Epelbaum E, Krebs H, Lee D and Meissner U G 2010 Phys.Rev.Lett. 104 142501 

Rosenstein B, Warr B and Park S 1991 Phys.Rept. 205 59-108 

Herbut I F, Juricic V and Roy B 2009 Phys. Rev. B 79 085116 

Armour W, Hands S and Strouthos C 2011 Phys.Rev. B84 075123 

Grabowska D, Kaplan D B and Nicholson A N 20 1 2 (Preprint\ 1208.5760) 

Wolff U 2008 Nucl.Phys. B789 258-276 (Preprint Wioi . 2872| 

Wolff U 2009 Nucl.Phys. B814 549-572 (Prepri nt\0812 . 0617$ 

Wenger U 2009 Phys.Rev. D80 071503 (Preprint \0812 . 3565) 

Scalettar R T, Scalapino D J, Sugar R L and Toussaint D 1987 Phys. Rev. B 36 8632-8641 
Duane S, Kennedy A D, Pendleton B J and Roweth D 1987 Phys. Lett. B195 216-222 
Clark M and Kennedy A 2007 Phys.Rev.Lett. 98 051601 (Preprint \hep- lat / 0608015| > 

Assad F F and Evertz H G 2008 Computational Many Particle Physics vol 739 (Lect. Notes Phys. Springer) chap 10, pp 
277-356 

Hands S, Kocic A and Kogut J B 1993 Annals of Physics 224 29 - 89 ISSN 0003-4916 URL |http : //wwTT] 
Isciencedirect ■ com/ science/art icle/pii/SOO 034 91 6837 103 90| 

Chandrasekharan S 2010 Phys.Rev. D82 025007 (Preprint ^ 1 ■ 5 7 3 6| 

Chandrasekharan S and Li A 201 1 PoS LATTICE2011 058 (Preprint \ 1 1 1 1 . 5 2 7~6) 

Chandrasekharan S and Li A 2012 Phys. Rev. Lett. 108(14) 140404 (Preprint \llll . 72 04) URL |http : / / 1 ink . aps . | 
|org/doi/10 . 1103/PhysRevLett ■ 108 . 140404| 

Chandrasekharan S 2012 Phys.Rev. D86 021701 (Preprint \l2 5 . 0084) 

Chandrasekharan S and Li A 2012 Phys.Rev. D85 091502 (Preprint \1202 . 6572) 

Chandrasekharan S and Wiese U J 1999 Phys. Rev. Lett. 83 3116-3119 

Troyer M and Wiese U J 2005 Phys.Rev.Lett. 94 170201 (Preprint \cond-mat / 0408370) 

Sharatchandra H S, Thun H J and Weisz P 1981 Nucl. Phys. B192 205 

Kogut J B and Susskind L 1975 Phys.Rev. Dll 395 

Banks T, Susskind L and Kogut J B 1976 Phys.Rev. D13 1043 



Hands S, Kim S and Kogut J B 1995 Nucl. Phys. B442 364-390 (Preprint hep-lat/9501037 



Hands S, Lucini B and Morrison S 2002 Phys. Rev. D65 036004 (Preprint hep-lat / 010 9001 



Christofi S and Strouthos C 2007 JHEP 05 088 (Preprint [hep - lat/0612031) 

Karkkainen L, Lacaze R, Lacock P and Petersson B 1994 Nucl. Phys. B415 781-796 (Preprint ftiep- lat/9310020) 
Hands S 1995 Phys. Rev. D 51 5816-5826 

Chandrasekharan S and Li A 2012 Phys. Rev. Lett. 108(14) 140404 URL jhttp : //link . aps . org/doi/10 .110371 

|PhysRevLett ■ 108 ■ 140404] 
Chandrasekharan S 2008 PoS LATTICE2008 003 
Chandrasekharan S and Wiese U J 201 1 (Preprint \Tl 08.2461) 



